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Fiber-reinforced  polymer  matrix  composite  materials  display  quite  complex  deformation  and  failure 
behavior  under  ballistic/blast  impact  loading  conditions.  This  complexity  is  generally  attributed  to  a 
number  of  factors  such  as  (a)  hierarchical/multi-length  scale  architecture  of  the  material  microstructure;  (b) 
nonlinear,  rate-dependent  and  often  pressure-sensitive  mechanical  response;  and  (c)  the  interplay  of  various 
intrinsic  phenomena  and  processes  such  as  fiber  twisting,  interfiher  friction/sliding,  etc.  Material  models 
currently  employed  in  the  computational  engineering  analyses  of  ballistic/blast  impact  protective  structures 
made  of  this  type  of  material  do  not  generally  include  many  of  the  aforementioned  aspects  of  the  material 
dynamic  behavior.  Consequently,  discrepancies  are  often  observed  between  computational  predictions  and 
their  experimental  counterparts.  To  address  this  problem,  the  results  of  an  extensive  set  of  molecular-level 
computational  analyses  regarding  the  role  of  various  microstructural/morphological  defects  on  the  Kevlar® 
fiber  mechanical  properties  are  used  to  upgrade  one  of  the  existing  continuum-level  material  models  for 
fiber-reinforced  composites.  The  results  obtained  show  that  the  response  of  the  material  is  significantly 
affected  as  a  result  of  the  incorporation  of  microstructural  effects  both  under  quasi-static  simple  mechanical 
testing  condition  and  under  dynamic  ballistic-impact  conditions. 


Keywords  ballistics,  composites,  Kevlar,  material  models,  micro- 
structural  defects 


1.  Introduction 

The  subject  of  the  present  study  is  high  specific- strength  and 
high  specific-stiffness  polymeric  fiber-reinforced  (e.g.,  Kevlar®, 
Twaron®,  etc.)  polymer-matrix  composite  materials  and  struc¬ 
tures.  These  materials/structures  are  commonly  used  in  various 
protective  systems  whose  main  requirement  is  a  high  level  of 
mass-efficient  penetration  resistance  against  high  kinetic  energy 
projectiles  (e.g.,  bullets,  mine,  lED  or  turbine  fragments,  etc.). 
Development  of  the  aforementioned  protection  systems  is 
traditionally  carried  out  using  legacy  knowledge  and  extensive 
fabricate-and-test  procedures.  Since  this  approach  is  not  only 
economically  unattractive,  but  is  often  associated  with 
significantly  longer  lead  times,  it  has  gradually  become 
complemented  by  the  appropriate  cost-  and  time-efficient 
computer-aided  engineering  analyses.  This  trend  has  been 
accelerated  by  the  recent  developments  in  the  numerical 
modeling  of  transient  nonlinear  dynamics  phenomena  such  as 
those  accompanying  ballistic  and  blast-loading  conditions. 
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However,  the  tools  used  in  these  analyses  themselves  suffer 
from  a  number  of  deficiencies/limitations  which  prevent  these 
analyses  from  being  more  widely  utilized.  One  of  the  main 
deficiencies  stems  from  the  inability  of  currently  available 
material  models  to  realistically  represent  the  response  of  these 
materials  under  high-deformation  rate,  large-strain,  and  high- 
pressure  loading  conditions,  the  conditions  typically  encoun¬ 
tered  during  projectile  impact  events.  One  of  the  reasons  for  the 
indicated  shortcomings  of  the  (continuum-type)  material  mod¬ 
els  currently  in  use  is  a  lack  of  inclusion  of  the  contributions  of 
various  phenomena  and  processes  occurring  at  different  length 
scales  to  the  overall  behavior  of  the  material  under  dynamic 
loading  conditions.  In  the  present  study,  an  attempt  is  made  to 
upgrade  an  existing  high-rate  composite  material  model 
through  incorporation  of  the  fine-scale  microstructural  (includ¬ 
ing  defect)  effects.  Hence,  the  two  main  topics  covered  in  the 
present  study  include  (a)  continuum-level  composite  material 
models;  and  (b)  multi-length  scale  character  of  the  composite- 
material  microstructure.  These  topics  are  briefiy  reviewed  in  the 
remainder  of  this  section. 

7.7  Continuum-Level  Composite  Material  Models 

A  review  of  the  open  literature,  carried  out  as  part  of  the 
present  study,  identified  number  of  continuum-level  material 
models  for  the  composite  materials  under  investigation.  Some 
of  these  material  models  (e.g..  Ref  1)  focus  on  the  role  of 
various  fiber/matrix  damage  processes  solely  due  to  in-plane 
loading.  Since  under  common  ballistic/blast-impact  loading 
conditions,  through-the-thickness  normal  and  transverse  shear 
stresses/strains  are  expected  to  play  an  important  role,  these 
types  of  models  are  of  little  utility  in  the  computational  analysis 
of  protective  structures.  Among  the  material  models  which 
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include  the  eontribution  of  all  the  attendant  stress/strain 
eomponents  on  the  material  response,  most  fall  into  the 
eategory  of  eontinuum  damage  meehanies  models.  In  other 
words,  damage  is  treated  as  a  eontinuum-filed  variable  whieh 
evolves  during  loading  and  degrades  loeally  the  material’s 
stiffness.  Among  these,  the  models  whieh  appear  most 
noteworthy  are  (a)  a  plane-stress  model  for  the  eomposite 
materials  reinforeed  with  unidireetionally  oriented  fibers 
(Ref  2);  (b)  a  fully  three-dimensional  eomposite-material  model 
whieh  aeeounts  for  the  progressive  eharaeter  of  the  main  fiber- 
damage  modes,  i.e.,  fiber  fraeture  under  axial  tension,  fiber 
buekling  under  eompressive  loading,  fiber  eutting  under 
transverse  shear,  and  fiber  erushing  under  eompressive 
through-the-thiekness  strains  (Ref  3,  4);  and  (e)  more  reeent 
three-dimensional  (3D)  eomposite  material  models  whieh 
inelude  various  eontinuously  distributed  statie  intralamina  and 
interlamina  damage  meehanisms  (e.g.,  fiber  breakage  within  the 
yams,  fiber/matrix  de-bonding,  diffuse  delamination/interlam¬ 
ina  separation,  ete.)  as  well  as  diserete  damage  modes 
(e.g.,  transverse  miero-eraeking)  (Ref  5-7). 

The  eomposite-material  models  briefiy  overviewed  above  all 
suffer  from  a  eommon  defieieney,  i.e.,  they  fail  to  aeeount  for 
the  effeets  of  high-loading  rates  under  ballistie/blast  impaet 
eonditions.  This  defieieney  is  removed  in  the  eomposite 
material  model  proposed  by  Yen  (Ref  8)  when  the  eontribution 
of  the  entire  strain/stress  field  to  all  the  main  fiber  and  matrix 
failure  modes  is  taken  into  aeeount.  This  model  is  adopted  in 
the  present  study,  and  an  attempt  is  made  to  upgrade  this  model 
to  inelude  the  effeets  of  fine  seale  mierostmeture  to  the  material 
response  under  dynamie  loading  eonditions.  Sinee  this  model  is 
overviewed  in  great  detail  in  seetion  2,  it  will  not  be  diseussed 
any  further  here. 

7.2  Hierarchical/Multi-Length  Scale  Material  Microstructure 

Filament-based  materials/stmetures  under  investigation  are 
quite  eomplex,  and  this  eomplexity  ean  be  linked  to  the 
following  main  sourees:  (a)  they  eontain  a  hierarehieal/multi- 
length  seale  arehiteeture;  (b)  their  meehanieal  response  is  often 
quite  nonlinear  and  rate/time-dependent;  and  (e)  they  are 
assoeiated  with  a  multiplieity  of  eomplex  phenomena/proeesses 
(e.g.,  filament  twisting,  interfilament  frietion  and  sliding,  ete.). 
In  our  reeent  studies  (Ref  9,  10),  an  attempt  was  made  to  help 
us  elarify  the  nature  of  the  multi-length  seale  hierarehy  of  the 
fabrie  struetures  or  fiber-reinforeed  eomposite  materials  under 
investigation.  A  set  of  eight  length  seales  (deseribed  later)  was 
established  and  briefiy  deseribed.  Sehematies  and  explanations 
of  the  set  of  eight  mierostruetural  length  seales  is  provided  in 
Fig.  1.  The  first  eolumn  in  this  figure  shows  a  set  of  simple 
sehematies  of  the  material  mierostrueture/arehiteeture  at  a  given 
length  seale  along  with  the  labels  used  to  denote  the  main 
mierostruetural  eonstituents.  In  the  seeond  eolumn,  a  brief 
deseription  is  provided  of  the  material  model  used  to  eapture 
the  material  behavior  at  the  length  seale  in  question.  The  key 
features  of  the  material  mierostrueture/arehiteeture  and  the 
eorresponding  material  model  at  eaeh  of  the  length  seales  listed 
in  Fig.  1  ean  be  summarized  as  follows:  (a)  Laminate  length- 
scale:  At  this  length  seale,  the  material  possesses  no  diseemable 
mierostruetural  features,  i.e.,  it  is  eompletely  homogenized.  An 
example  of  the  laminate-length  seale  eomposite  material  model 
ean  be  found  in  Ref  1 1 ;  (b)  Stacked-lamina  length-scale:  At  this 
length  seale,  the  presenee  of  diseretely  staeked  laminae  is 
reeognized,  while  the  material  within  eaeh  lamina  as  well  as 


interlamina  boundaries  are  kept  featureless/homogenized.  An 
example  of  the  eomposite  material  model  at  this  length  seale 
ean  be  found  in  Ref  7;  (e)  Single-lamina  length-scale:  In  this 
ease,  the  mierostrueture/arehiteeture  of  the  reinforeement  and 
matrix  phases  are  explieitly  taken  into  aeeount,  while  the  two 
assoeiated  materials  are  treated  as  featureless/homogenized. 
The  resulting  '‘two-phase''  eomposite  lamina  model  for  eaeh 
lamina  is  then  eombined  with  a  fully  homogenized  lamina/ 
lamina  interfaee  model  to  form  a  staekedlamina  eomposite 
laminate  strueture.  An  example  of  this  type  of  material  model 
ean  be  found  in  Ref  12;  (d)  Fabric  unit-cell  length-scale:  At 
this  length  seale,  a  eloser  look  is  given  to  the  arehiteeture  of  the 
woven  fabrie.  Speeifieally,  details  of  yam  weaving  and 
erimping,  yam  eross-seetion  ehange,  and  yam  sliding  at  the 
warp-yam/weft-yam  eross-over  points  are  taken  into  aeeount. 
Eaeh  yam  is  fully  homogenized  in  this  ease  as  are  the  matrix 
and  the  lamina/lamina  interfaees.  An  example  of  the  model  of 
this  type  ean  be  found  in  Ref  13-15;  (e)  Yarn  length-scale:  In 
this  ease,  the  internal  stmeture/arehiteeture  of  eaeh  yam  is 
aeeounted  for  explieitly.  In  other  words,  yams  are  eonsidered  as 
assemblies  of  nearly  parallel  fibers/filaments  whieh  are 
meehanieally  engaged  by  either  the  applieation  of  a  light  twist 
to  the  yam  or  by  wrapping  a  fiber  around  the  fiber/filament 
assembly.  Typieally,  the  detailed  mierostmeture  of  the  yams  is 
not  ineorporated  into  material  models  used  in  large-seale 
simulations  of  projeetile/proteetion- system  impaet  interaetions 
(due  to  the  unmanageably  large  eomputer  resourees  needed). 
Instead,  eomputational  results  of  the  meehanieal  response  of  the 
individual  yams,  when  subjeeted  to  a  variety  of  loading 
eonditions,  are  used  to  enrieh  the  homogenized  material  model 
for  the  yams.  The  latter  is  subsequently  used  at  the  fabrie  unit- 
eell  length  seale.  An  example  of  this  approaeh  ean  be  found  in 
Ref  16;  (f)  Fiber  length-scale:  At  this  length  seale,  fibers  are 
eonsidered  as  assemblies  of  aligned  long-ehain  moleeules 
whieh  are  held  together  by  non-bonding  (van  der  Waals  or 
Coulomb)  forees.  In  the  ease  of  Kevlar®  type  fibers,  there  is  a 
substantial  experimental  support  (e.g..  Ref  17)  for  the  existenee 
of  fibrils  within  the  fibers.  Fibrils  are  smaller  bundles  of 
moleeular  ehains  within  whieh  ehain  moleeules  are  tightly 
bonded  into  a  perfeet  or  nearly  perfeet  erystalline  phase.  Thus, 
in  the  ease  of  Kevlar®  type  fibers,  fibers  ean  be  eonsidered  as 
an  assembly  of  fibrils;  (g)  Fibril-level  length-scale:  As 
mentioned  above,  within  the  fibrils,  the  moleeular  stmeture  is 
erystalline  or  nearly  erystalline.  However,  as  elearly  demon¬ 
strated  in  our  reeent  study  (Ref  9),  the  material  at  this  length 
seale  (as  well  as  the  fiber  length  seale)  may  eontain  a  variety  of 
mierostmetural  and  topologieal  defeets  and  ehemieal  impurities 
whieh  may  signifieantly  alter  its  properties.  As  in  the  fiber 
length  seale  ease,  the  material  is  treated  as  a  eolleetion  of 
diserete  interaeting/bonded  partieles  and  analyzed  using  the 
atomie/moleeular  modeling  tools/proeedures;  and  (h)  Molecu¬ 
lar  chain-level  length-scale:  At  this  length  seale,  ehemieal 
stmeture  and  eonformation  of  the  individual  moleeules  eonsti- 
tuting  the  ehain  are  analyzed  using  the  aforementioned  atomie/ 
moleeular  modeling  tools/proeedures.  The  main  goal  of  the 
analysis  at  this  length  seale  is  to  identify  the  most-likely 
moleeular  eonformations  present  in  the  fibrils. 

7.3  Objective 

As  mentioned  earlier,  the  main  objeetive  of  the  present  study 
is  to  identify  and  quantify  fine-seale  mierostmetural  effeets  on 
the  dynamie  meehanieal  response  of  the  eomposite  materials 
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Fig.  1  Various  length  scales  and  the  associated  material  model  assumptions/simplifications  used  in  the  study  of  polymer-matrix  composite 
materials  with  high-performance  fiber-based  structures 


under  investigation  and  to  ineorporate  these  effeets  into  the 
eontinuum-level  material  model  of  Yen  (Ref  8). 

1.4  Organization  of  the  Paper 

A  brief  overview  of  the  eontinuum-level  eomposite  material 
model  proposed  by  Yen  (Ref  8)  is  provided  in  seetion  2.  A 
summary  of  the  fiber  and  sub-fiber  level  mierostrueture  and 
morphology  of  Kevlar® -based  materials  ineluding  the  polymer- 
synthesis  and  fiber  fabrieation  proeesses-indueed  defeets  is 
presented  in  seetion  3.  The  proeedure  used  to  upgrade  the 
model  of  Yen  (Ref  8)  and  the  results  obtained  using  the 
upgraded  material  model  are  all  presented  and  diseussed  in 
seetion  4.  The  key  findings  resulting  from  the  present  study  are 
summarized  in  seetion  5. 


2.  Overview  of  the  Yen  Composite  Materiai 
Modei  (Ref  8) 

In  this  seetion,  a  brief  overview  is  provided  of  the  rate- 
dependent  eomposite  material  model  developed  by  Yen  (Ref  8). 
The  defining  features  of  this  model  ean  be  summarized  as 
follows: 

(a)  The  model  is  an  extension  of  the  Hashin  eomposite- 
material  failure  model  (Ref  18)  whieh  more  aeeurately 
aeeounts  for  the  three-dimensional  nature  of  the  stress/ 
strain  states,  high  strain  rates,  and  high  pressures 
eneountered  during  ballistie  impaet  of  eomposite  lami¬ 
nates  reinforeed  with  plain-woven  fabrie; 


Journal  of  Materials  Engineering  and  Performance 


Volume  22(3)  March  2013—683 


Yarn  Length-scale 

Eich  yam  ts  consrdered  as  an  assembly 
of  f ully-homoge  nized  mecha  nica  I  ly* 
engaged  discrete  fibers 
Rber-fiber  contact  and  sliding  are 
modeled  explicitiy 


w*  ' 


Fiber  Length-scale 

Rbers  are  considered  as  an  assembly  of 
nearly-parallel  fibrils  held  together  by 
weak  van  der  Walls  forces 
Adjacent  fibrils  differ  with  respect  to  the 
orientation  of  the  stacked  sheets. 


Fibril  Length-scale 

Rbrils  are  considered  as  crystalline 
thread-like  entities  containing  parallel 
sheets  of  aligned  and  hydrogen 
bonded  molecular  chains, 
inter-sheet  bonding  is  of  a  van  der  Walls 
character. 


Molecular-level  Length-scale 

«  Individual  molecules  are  represented  as 
assemblies  of  covalently  bonded 
atoms. 

*  The  effect  of  hydrogen  bonding  on  the 
molecular  topology  (e  g,  the  formation 
of  hetical  structure  is  considered 
explicitly). 


Fig.  1  continued 


(b)  The  material  model  is  strain  based  and  utilizes  the  eon- 
eept  of  a  ply/lamina-averaged  loeal  strain; 

(e)  The  model  ineorporates  Mohr-Coulomb  (Ref  19)  shear 
strength  effeets,  i.e.,  it  aeeounts  for  the  faet  that  the 
(matrix)  in-plane  shear  strength  is  generally  a  funetion 
of  the  (tensile/eompressive)  eharaeter  and  magnitude  of 
the  through-the-thiekness  lamina  strain; 

(d)  Within  the  model,  no  distinetion  is  made  between  fibers 
and  yams.  Henee,  to  be  eonsistent  with  the  original 
model,  in  this  seetion,  the  term  '‘fiber''  will  be  used  to 
denote  the  thread-like  entity  whieh  is  used  in  the  fabrie- 
weaving  proeess.  In  the  next  and  the  subsequent  see- 
tions,  the  same  entity  will  be  referred  to  as  a  "yarn"; 

(e)  The  following  three  fiber  and  two  matrix  failure  modes 
are  eonsidered:  (i)  fiber  failure  due  to  a  eombination  of 
axial  tensile  and  transverse  shear  strains;  (ii)  fiber  buek- 
ling  failure  due  to  axial  eompressive  strains;  (iii)  fiber 
emshing  failure  due  to  lamina  through-the-thiekness 
eompressive  strains;  (iv)  matrix  in-plane  shear  failure; 


and  (v)  matrix  delamination  failure  due  to  a  eombined 
effeet  of  the  through-the-thiekness  tensile  strains  and 
transverse  shear  strains;  and 

(f)  The  model  is  essentially  of  a  linear-elastie  orthotropie 
type  with  progressively  degradable  elastie  stiffness  eon- 
stants.  In  other  words,  a  generalized  Hooke’s  law  relation 
is  assumed  to  hold  between  stress  and  strain  quantities, 
while  damage  is  treated  as  a  proeess  whieh  takes  plaee 
gradually  during  ballistie  impaet.  Consequently,  the  mod¬ 
el  relies  on  the  definition  of  the  damage  initiation  eriteria 
and  the  eorresponding  damage  evolution  laws. 

2.7  Damage  Initiation  Criteria 

In  this  seetion,  a  brief  overview  is  provided  of  the  funetional 
relations  defining  the  onset  (and  eontinuation)  of  damage  for 
eaeh  of  the  five  aforementioned  fiber/matrix  failure  modes.  For 
eaeh  of  the  failure  modes  in  question,  as  will  be  shown  shortly, 
beeause  distinetion  is  made  between  weft  and  warp  fibers,  a  set 
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of  seven  equations  had  to  be  used  to  fully  deseribe  damage 
initiation.  Eaeh  of  the  seven  equations  has  a  general  form: 
f  —  rf  =  0  7)  where  /  represents  the  appropriate 

strain-based  loading  funetion,  while  represents  the  eorre- 
sponding  strain-based  damage  threshold.  Initially,  eaeh  is  set 
to  a  value  of  1.0  and,  as  a  result  of  progressive  damage, 
values  eontinue  to  inerease. 

It  should  be  noted  that  in  the  remainder  of  this  artiele,  the 
following  Cartesian  eoordinate  system  is  used:  x  and  y  are, 
respeetively,  aligned  with  the  weft  and  the  warp  fibers,  while  z  is 
aligned  in  the  eomposite- laminate  through-the-thiekness  direetion. 

2.1.1  Fiber  Damage  Under  Combined  Axial  Tension  and 
Transverse  Shear.  The  onset  and  progression  of  this  failure 
mode  is  assumed  to  be  defined  by  a  superposition  of  the 
quadratie  eontributions  of  the  axial  tensile  and  transverse  shear 
strains.  For  the  weft  and  the  warp  fiber  direetions,  the 
eorresponding  damage  initiation  funetions  ean  be  written  as 

^  (^)  =0  (Eq  1) 


fi-fi 


I  fGyzEyz\  _  n 


(Eq  2) 


where  E  denotes  the  fiber  Young’s  modulus,  G  the  eorre¬ 
sponding  shear  modulus,  S  the  fiber  strength,  s  the  relevant 
strain  eomponent,  the  Maeaulay  braekets  ()  are  used  to  de¬ 
note  the  ply/lamina-average  loeal  quantity,  and  subseripts  T 
and  FS  denote,  respeetively,  tension  and  fiber  shear.  It  should 
be  noted  that  in  Fq  1  and  2,  the  laminate  transverse  shear 
strength  is  assumed  to  be  eontrolled  by  the  eorresponding 
shear  strength  of  the  fiber  itself 

2.1.2  Fiber  Damage  Under  Axial  Compression.  The 
onset  and  progression  of  this  failure  mode  is  assumed  to  be 
defined  by  the  eompressive  axial  strains,  effeet  of  whieh  is 
modulated  by  the  lamina  through-the-thiekness  eompressive 
strain.  In  other  words,  axial  eompression-indueed  fiber  buek- 
ling  is  resisted  by  the  stabilizing  effeet  of  the  through-thiekness 
eompressive  strain.  For  the  weft  and  the  warp  fiber  direetions, 
the  eorresponding  damage  initiation  funetions  ean  be  written  as 


h-ri  = 


SxC 


—rj  =  0 


(Fq3) 


f4->i  = 


(Eq  4) 


where  prime  is  used  to  denote  modulation  of  the  axial  eom¬ 
pressive  strains  as  =  —  S;^  —  {—^z){Ez/Ex)  and  =  ~^y~ 
{—£.z){Ez/Ey),  and  the  subseript  C  denotes  eompression. 

2.1.3  Fiber  Damage  Under  Transverse  Compressive 
Loading.  This  failure  mode  involves  fiber  erushing  under 
sufheiently  high  through-the-thiekness  eompressive  strain,  and 
its  onset/progression  is  governed  by  the  following  equation: 

=  0  (Eq5) 

where  the  subseript  FC  denotes  a  fiber  eompression  quantity. 

2.1.4  Matrix  In-plane  Shear  Failure.  This  mode  of 
failure  is  assoeiated  with  the  faet  that  the  matrix  shear  strength 
is  lower  than  the  fiber  shear  strength  and  that  the  eomposite 


material  in  question  possesses  a  laminated  arehiteeture.  Con¬ 
sequently,  sufheiently  high  in-plane  shear  strains  may  eause 
matrix  shear  failure  without  giving  rise  to  any  signifieant  fiber 
damage.  The  onset  and  progression  of  this  type  of  damage  are 
governed  by  the  following  funetional  relationship: 

/6-^=(^)'-^  =  0  (Eq6) 


2.1.5  Matrix  Delamination  due  to  Through-the-Thick- 
ness  Tensile  and  Shear  Stresses.  This  mode  of  failure  also 
gives  rise  to  the  matrix  in-plane  failure,  and  its  onset/ 
progression  is  deseribed  by  the  following  relation: 


,  SzT  )  ySyzoESsR/  \SzxOESsr)  j 


=  0 


(Eq  7) 


Examination  of  Eq  7  reveals  the  presenee  of  three  terms  within 
the  eurly  braekets  on  its  right-hand  side.  The  first  term  defines  the 
eontribution  of  tensile  through-the-thiekness  strains,  while  the 
remaining  terms  represent  the  eontributions  of  the  transverse  shear 
strains,  to  matrix  delamination.  In  Eq  7,  subseript  0  is  used  to 
denote  the  appropriate  shear  strength  in  the  absenee  of  the 
assoeiated  eompressive  normal  strain,  while  subseript  SR  is  used  to 
denote  an  inerease  in  this  shear  strength  due  to  the  presenee  of  a 
normal  eompressive  strain.  The  ^Ssr  (sliding  resistanee  term  is 
assessed  using  the  Mohr-Coulomb  approaeh  as 

•SsR  =  tan  (|)(-Ez)  (Eq  8) 

where  tan(|)  is  the  frietion  angle. 

Equation  7  also  eontains  a  stress-eoneentration  (sealing) 
faetor,  S.  This  faetor  is  assigned  a  value  larger  than  1.0  in  the 
non-delaminated  portion  of  the  material  adjaeent  to  the 
delamination  boundary/front.  This  algorithm  has  been  found 
to  introduee  enough  flexibility  into  the  model  so  that  a  better 
agreement  ean  be  obtained  with  the  eorresponding  experimen¬ 
tal  results  pertaining  to  the  temporal  evolution  of  the  delam¬ 
ination  zone  (Ref  8); 


2.2  Damage  Evolution  and  Material  Degradation 

2.2.1  Damage  Surface.  Equations  1-7  given  above 
define  within  the  strain  spaee  a  damage  initiation/progression 
surfaee.  Examination  of  these  equations  reveals  that  the 
material  stiffness  is  eharaeterized  by  three  Young’s  moduli 
and  three  shear  moduli.  The  initial  damage  initiation/progres¬ 
sion  surfaee  is  defined  by  substituting  the  initial  values  of  the 
six  moduli  into  Eq  1-7.  As  the  material  moduli  degrade,  the 
failure  surfaee  within  the  strain  spaee  expands  (i.e.,  additional 
straining  is  required  to  produee  additional  damage).  Within  the 
eomposite  material  model  developed  by  Yen  (Ref  8),  damage 
progression  is  assumed  to  eause  degradation  in  material 
stiffness  (without  an  aeeompanying  loss  in  material  strength) 
properties.  The  same  assumption  will  be  adopted  in  the  present 
study.  In  our  fortheoming  study,  however,  it  is  planned  to 
address  the  problem  of  strength  degradation.  It  should  be  noted 
that  although  stiffiiess  is  eonsidered  as  a  non-degradable 
material  property,  it  will  be  eonsidered  as  a  stoehastie  rather 
than  a  deterministie  quantity  in  the  present  study. 

2.2.2  Stiffness  Degradation  Variables.  To  aeeount  for 
the  aforementioned  damage-indueed  material- stiffness  degra¬ 
dation,  six  stiffness  degradation  variables,  coj  {j  =  1,...,6),  are 
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first  introduced.  Then,  the  effeet  of  damage  on  the  material 
elastie  eomplianee  matrix  is  deseribed  by  the  following  relation 
(Ref  2): 


[5]  = 


1  -y 

(1-(Di)£'jc  Ey 

-y  1 

{\-(£>2)Ey 
-Vxz  -Vjz 

Ex  Ey 

0  0 

0  0 

0  0 


-Vzx 

E, 

E, 

1 

(l-yR 


0 

0 
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0  0 

0  0 

0  0 

(1-CD4)G^  ^ 

0  _ 1 _ 

0  0 


0 

0 

0 

0 

0 

1 

(1-(D6)Gzjc 


(Eq9) 


Finally,  the  ''degraded''  material  elastie  stiffness  matrix  is 
obtained  by  inverting  the  material  elastie  eomplianee  matrix  as 
defined  by  Eq  9.  This  proeedure  shows  that  to  eompute  the 
instantaneous  (degraded)  material  stiffness  properties,  one  must 
know  the  eorresponding  values  of  the  six  stiffness  degradation 
variables.  To  evaluate  the  same,  one  must  define  the  appropriate 
evolution  equations  for  the  stiffness  degradation  variables. 

Following  Matzenmiller,  et  al.  (Ref  2),  the  needed  evolution 
equations  are  defined  as 

=  (Eq  10) 

i 

where  (j)^  (/  =  l,.--?7)  denotes  the  rate  of  progression  of  eaeh 
of  the  failure  modes  (the  initiation  of  whieh  was  deseribed  by 
Eq  1-7),  while  qji  represents  a  binary  (0  or  1)  two  dimen¬ 
sional  array  with  j  =  1,...,6  and  quantifies  the  eontribution  of 
different  failure  modes  to  the  evolution  of  stiffness  degrada¬ 
tion  variables. 

2.2.3  Damage  Evolution  Laws.  Next,  (|)^  fimetional  rela¬ 
tions  are  derived  for  eaeh  of  the  failure  modes.  It  should  be  first 
noted  that  failure  and  the  assoeiated  stiffness  degradation  are 
irreversible  and  that,  during  loading,  an  inerease  in  the  extent  of 
damage  of  a  given  mode  oeeurs  only  if  the  resulting  strain 
inerement  erosses  the  eorresponding  failure  surfaee,  i.e.,  that 
the  projeetion  of  the  strain  inerement  onto  the  loeal  failure 
surfaee  outward  normal  is  positive  as  Y^k{dfi/dEk)Ek>  0, 
where  k  (=  1,...,6)  is  used  to  denote  six  eomponents  of  the 
strain  veetor.  Then,  by  defining  the  failure-mode-speeifie 
progression  funetion,  the  rate  of  progression  of  failure 
mode  i  ean  be  defined  as 


(Eq  11) 


where  no  summation  is  implied  over  repeating  indiees.  By 
ehoosing 

Y,=^(l-rf’'  (Eql2) 

where  the  material  parameter  m  quantifies  the  stiffness  sensi¬ 
tivity  to  damage  extent  and  noting  that  {dfi/dSk)Ek  —  f 
the  final  expression  for  the  rate  of  failure  progression  ean  be 
defined  as 


(j)^- =  -  (1  —  (j)^)^^  ^fi  (no  summation  over /)  (Eq  13) 

2.2.4  Derivation  of  the  Damage  Coupling  Matrix.  To 

assemble  the  damage  eoupling  matrix,  qj^  (j  denotes  the  elastie 
moduli,  while  i  denotes  the  failure  mode),  the  following 
approaeh  was  taken: 


(a)  Failure  of  the  weft  and  warp  fibers  both  in  tension  and 
transverse  shear  and  in  eompression  degrades  the  assoei¬ 
ated  Young’s  modulus  and  shear  moduli; 

(b)  Fiber  erushing  mode  degrades  all  six  elastie  moduli; 

(e)  In  the  ease  of  the  in-plane  shear  failure,  only  the  assoei¬ 
ated  shear  modulus  is  affeeted;  and 
(d)  Delamination  degrades  the  through- the-thiekness  Young’s 
modulus  as  well  as  the  two  transverse  shear  moduli. 


The  resulting  damage  eoupling  matrix  is 


0  1  0 
1  0  1 
0  0  0 
1  1  1 
1  0  1 
0  1  0 


1  0  0 
1  0  0 
1  0  1 
1  1  0 
1  0  1 
1  0  1 


(Eq  14) 


2.2.5  Integrated  Form  of  the  Stiffness  Degradation 
Variables.  By  eombining  Eq  10,  12,  and  13,  the  following 
expression  is  obtained  for  the  eontribution  of  the  failure  mode, 
z,  to  eaeh  stiffness  degradation  variable: 

CO,.  =  1 -e»(' ■'!”),  ri>\  (Eq  15) 

Examination  of  Eq  15  reveals  that  initially,  when  ri  is  equal 
to  1.0,  the  eorresponding  eontribution  to  cuy  is  zero.  During 
loading,  when  the  strain  inerement  eauses  the  strain  to  fall 
outside  the  damage  surfaee  assoeiated  with  failure  mode  z,  the 
eorresponding  damage  surfaee  erossed  is  expanded  outward. 
This,  in  aeeordanee  with  Eq  1-7,  gives  rise  to  an  inerease  in 
and,  in  turn,  via  Eq  15  to  an  inerease  in  cuy. 

2.2.6  Residual  Stiffness.  Aeeording  to  Eq  15,  at  sufh- 
eiently  high  r^,  the  stiffness  degradation  variable  may  aequire  a 
value  of  1.0  in  whieh  ease  the  material  eompletely  loses  the 
assoeiated  load  earrying  eapaeity.  While  this  is  physieally 
realistie  for  some  of  the  failure  modes,  e.g.,  the  fiber  failure 
under  axial  tension/transverse  shear,  it  may  not  be  fully  realistie 
under  other  failure  modes,  e.g.,  fiber  axial  eompression. 
Consequently,  the  following  modifieations  are  made  in  Eq  15: 


(a)  For  the  fiber  axial  eompression  failure  mode: 
(Dy  =  ^1  —  (1  —  r|^),  z  =  3  or4 


(Eq  16) 


where  the  subseript  R  is  used  to  denote  the  respeetive  resid¬ 
ual-strength  quantities;  and 


(b)  In  the  ease  of  delamination  failure  mode,  distinetion  is 
made  between  this  failure  mode  under  eompressive  and 
tensile  through- the-thiekness  strains.  In  the  former  ease, 
delamination  surfaees  are  assumed  to  be  elosed,  and 
Eq  15  for  the  two  transverse  shear  moduli  is  modified 
to  inelude  the  effeet  of  the  sliding  resistanee  as 

^SR 

(Eq  17) 
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.  On  the  other  hand,  in  the  ease  of  the  delamination  under 
tensile  through-the-thiekness  strains,  no  modifieations  are  made 
to  Eq  15. 


2.3  Rate  Dependency  of  Strength  and  Stiffness 


As  mentioned  earlier,  one  of  the  defining  features  of  the 
model  developed  by  Yen  (Ref  8)  is  the  inelusion  of  material 
(strength  and  stiffness)  rate  dependeney.  To  aeeount  for  rate 
dependeney  of  different  eomponents  of  material  strength,  the 
following  sealing  relation  was  used: 

{Skt}  =  {^o}  (^1  +  Cl  (Eq  18) 

where 


✓ 

SyT 

SxC 

<  SyC 

{g}  =  < 

Sfc 

Szl 

SxFS 

^zx 

^yz 

Cl  is  the  strength  rate  eonstant,  and  the  6  x  1  eolumn  veetor 
{Sq}  is  equal  to  {-S'rt}  when  eaeh  of  the  relevant  strain  rates 
is  equal  to  the  single  referenee  strain  rate  Sq  =  1.0  s“^. 

The  rate  dependeney  of  different  stiffness  eomponents  was 
aeeounted  for  in  a  similar  fashion  as 

{£„}  =  {£o}  (^1  +  C2  (Eq  19) 

where 


{^rt}  = 


Ey 

E, 


>, 


C2  is  the  stiffness  rate  eonstant,  and  the  6  x  1  eolumn  veetor 
{E’o}  is  equal  to  {E’rt}  when  eaeh  of  the  relevant  strain  rates 
is  equal  to  the  single  referenee  strain  rate  Sq  =  1.0  s“^. 


2.4  Material  Model  Parameterization 

The  material  model  presented  above  eontains  14  strength 
parameters  Y^fs^  ^yFS>  ^xc>  ^xy  ^zv  ^yzO>  ^zxO> 

phi,  YcCr?  and  SyCR),  six  stiffness  parameters  (E^,  Ey,  E^, 

Gyz,  and  Gz^\  two  rate  dependeney  parameters  (Ci  and  C2),  and 
two  damage  parameters  {S  and  m).  These  parameters  are 
eommonly  assessed  by  applying  regression  analysis  to  the 
results  from  a  variety  of  meehanieal  tests  over  a  range  of  strain 
rates.  The  large  number  of  the  model  parameters  makes  this 
parameterization  effort  quite  eumbersome,  time  eonsuming  and 
eostly.  In  addition,  as  will  be  diseussed  in  the  next  seetion, 
some  of  these  parameters  are  stoehastie  quantities  and  their 
eomplete  speeifieation  requires  the  use  of  the  appropriate 
probability  distribution  funetions.  Determination  of  these 
distribution  funetions  is  even  more  eostly.  Henee,  a  multi¬ 
length  seale  eomputational  approaeh  is  introdueed  in  the 
next  seetion  whieh  allows  establishment  of  the  appropriate 


distribution  funetions  through  the  use  of  time-efheient,  low  eost 
eomputational  methods  and  tools. 


3.  Microstructure/Property  Relations  in  Kevlar® 
Fibers 

In  this  seetion,  a  brief  overview  is  provided  of  the  basie 
mierostruetural  features  in  Kevlar®  at  the  fibril/fiber/yam  length 
seale(s),  and  their  effeets  on  the  key  material  properties. 

3.1  Fiber  Crystalline  Structure 

Kevlar®  (p-phenylene  terephthalamide,  PPTA)  belongs  to 
the  family  of  polymerie  materials  known  as  polyamides. 
Polyamides  are  typieally  elassified  as  aromatie  polyamides  or 
aramids  (e.g.,  Kevlar®,  the  subjeet  of  the  eurrent  study)  and 
non-aromatie  polyamides  (e.g..  Nylon  6,6).  Kevlar®  ehains  are 
inherently  stiff  and  henee  the  ehains  do  not  readily  flex. 
Consequently,  in  eontrast  to  the  flexible  polymerie  moleeules 
whieh  ean  undergo  extensive  folding  and  form  the  eommon 
(erystalline  +  amorphous)  two  phase  strueture,  the  Kevlar® 
flbers  tend  to  typieally  form  either  a  paraerystalline  or  a  fully 
erystalline  strueture.  The  tendeney  for  the  formation  of 
paraerystalline  or  erystalline  struetures  is  promoted  by  the 
presenee  of  the  planar  phenylene  and  amide  groups  and  by 
the  ability  of  the  adjaeent  ehains  to  form  hydrogen  bonds.  In  the 
ease  of  the  paraerystalline  strueture,  Kevlar®  moleeules  are  all 
aligned  in  the  same  direetion  but  no  order  exists  in  a  plane 
orthogonal  to  this  direetion.  In  sharp  eontrast,  in  the  ease  of  the 
fully  erystalline  Kevlar®  flbers,  moleeules  are  aligned  in  all 
three  mutually  orthogonal  direetions. 

The  Kevlar®  erystal  strueture  is  layered  and  eonsists  of 
parallel  (ABABAB...)  staeked  sheets  formed  due  to  the 
hydrogen  bonding  between  the  adjaeent  parallel  Kevlar® 
moleeules.  The  intersheet  bonding  on  the  other  hand  is  mainly 
of  the  van  der  Waals  (and  p-eleetron  weak  ehemieal-bond)  type. 
Owing  to  the  relatively  low  strength  of  the  van  der  Waals  and 
the  p-eleetron  type  weak  intersheet  bonding,  Kevlar®  flbers  are 
prone  to  the  formation  of  staeking  faults  and  kink  bands,  and 
typieally  possess  an  inferior  eompressive  strength  and  buekling 
resistanee. 

3.2  Common  Defects  in  Kevlai^  Fibers 

As  in  most  engineering  materials,  properties  of  Kevlar® 
flbers  are  greatly  affeeted  by  the  presenee  of  various  defeets/ 
flaws.  These  flaws  ean  be  introdueed  during  Kevlar®  synthesis 
and  Kevlar®  flber  fabrieation  proeesses.  In  our  reeent  study 
(Ref  10),  a  detailed  overview  was  provided  of  the  main  steps 
assoeiated  with  Kevlar®  synthesis  and  Kevlar®  flber  proeess- 
ing.  In  the  next  seetion,  a  brief  summary  is  provided  of  the 
eommonly  observed  synthesis/proeessing-indueed  defeets  in 
Kevlar®  fibers. 

A  summary  of  the  Kevlar®  fiber’s  most  eommon  defeets, 
their  dimensionality,  their  eause,  methods  of  redueing  their 
density,  and  their  typieal  eoneentrations  is  provided  in  Table  1 . 
Moleeular-level  sehematies  of  these  defeets  are  displayed  in 
Fig.  2(a)-(h).  These  defeets  are  expeeted  to  have  a  profound 
effeet  on  the  fiber  properties,  as  well  as  on  the  properties  of 
eoarser-seale  materials  and  struetures  eontaining  Kevlar® 
fibers  (e.g.,  yams,  fabries,  plies,  lamina,  and  laminates). 
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Table  1  A  summary  of  the  most  common  defects  found  in  Kevlar®-based  materials 


Flaw  group 

Flaw 

Causes 

How  causes  can  be  remedied/ 
mitigated 

Range  for  concentration 

Isolated  chain 

-COOH 

H2SO4  eatalyzed  hydrolysis 

Use  eoneentrated  H2SO4  for 

0.35  per  PPTA  ehain 

ends  (point 

eausing  PPTA  ehain  scission. 

dope  preparation.  Shorten  the 

for  eaeh  flaw(a) 

defeet) 

Na^  deficiency  with  respect 
to  complete  neutralization  of 
side/end  aeidie  groups 

fiber  wash  time 

(~  350  ppm-mass- 
based) 

-NH2 

H2SO4  eatalyzed  hydrolysis 
eausing  PPTA  ehain  seission. 
Na^  defieieney  with  respeet 
to  eomplete  neutralization  of 
side/end  aeidie  groups 

Use  higher  eoneentration 

NaOH  solution 

0.35  per  PPTA  ehain 
for  eaeh  flaw(a) 

(~  350  ppm-mass- 
based) 

-COO“Na+ 

COOH  neutralization  with  Na^ 

No  remedy  required  sinee 
this  is  one  of  the  preferred 
ehain  ends 

1 . 1  per  PPTA  ehain(a) 

(~  1100  ppm-mass- 
based) 

-NH3^HS04“ 

Sulfonation  of  the  NH2  ehain 
ends 

Inerease  the  H2SO4  removal 
and  neutralization  rate 

0.2  per  PPTA  ehain(a) 

(~  200  ppm-mass- 
based) 

Side  groups 

-SO3H 

Exposure  of  PPTA  in  the  dope 

Reduee  the  H2SO4  eoneen¬ 

~1300  ppm 

(point 

defeet) 

to  eoneentrated  H2SO4 
(sulfonation) 

tration  in  the  dope 

(mass-based) 

-S03“Na^ 

Neutralization  of  sulfonic  acid 
side  groups  by  NaOH 

Remedy  may  not  be  required 
sinee  this  side  group 
improves  fiber  longevity. 
However  meehanieal  perfor- 
manee  may  be  eompromised 

~2500  ppm 
(mass-based) 

Voids  and 

Mierovoids 

Swelling  indueed  by  hydration 

Inerease  the  extent  of  sodium 

~150  ppm 

interstitials 
(point  defeets) 

of  intra-fibrillar  Na2S04 

salt  dissolution  by  prolonged 
exposure  of  fibers  to  boiling 
water 

(mass-based) 

Mobile  trapped  H2SO4 

Non-neutralized  or  unwashed 
intra-fibrillar  H2SO4 

Thorough  washing  in  hot 
solvent  aqueous  bath 

~70  ppm  (mass-based) 

Defeet  bands 

NH3^HS04“ 

Coulombie  forees-indueed 

The  phenomenon  is  not 

One  band  every 

(planar  defeets) 

agglomerated 

elustering  of  ion-terminated 

well  understood  so  no 

40-60  nm  of  fibril 

(a)  Extruded  fibers 

ehain  ends 

ehain  ends 

remedy  is  obvious 

(ea.  3000  ppm- 
mass-based) 

Consequently,  a  high-fidelity  material  model  for  Kevlar®  fiber- 
reinforeed  eomposite  materials  must  aeeount  for  the  presenee, 
distribution,  and  poteney  of  these  defeets  in  modifying  material 
response  under  dynamie  loading  eonditions.  As  mentioned 
earlier,  the  main  objeetive  of  the  present  study  is  to  upgrade  one 
of  the  eurrently  used  eontinuum-level  eomposite  material 
models  via  ineorporation  of  the  effeet  of  these  defeets. 

3.3  Effect  of  Defects  on  Kevlai^  Properties 

In  this  seetion,  a  brief  summary  is  provided  of  the  main 
results  obtained  in  our  reeent  study  (Ref  10)  in  whieh 
moleeular-level  modeling  was  used  to  establish  the  effeets  of 
various  defeets  on  the  material  properties  (primarily  strength 
and  stiffness).  In  Ref  10,  a  proeedure  was  introdueed  to 
establish  relationships  between  topology/kinematies  and  prop¬ 
erties  of  Kevlar®  fibrils,  fibers,  and  yams.  Fig.  3(a)  and  (b).  In 
Fig.  3(a),  a  Kevlar®  fiber  is  depieted,  and  it  is  shown  that  it 
eonsists  of  a  number  of  nearly  parallel  fibrils.  Within  eaeh  of 
the  fibrils,  the  fibril’s  axis  lies  within  the  Kevlar®  sheets,  and 
the  fibrils  merely  differ  in  the  orientation  of  the  sheet  normals 
(the  sheet  normal  in  eaeh  fibril  is  orthogonal  to  the  fibril  axis). 
Figure  3(b)  displays  the  assumed  relationship  between  the 
fibers  and  the  yams.  Simply  stated,  the  yams  are  treated  as 


assemblies  of  nearly  parallel  fibers  whieh  are  given  a  light  twist 
to  meehanieally  engage  the  fibers. 

Within  the  analysis  presented  in  Ref  10,  it  was  assumed  that 
the  longitudinal  stiffness,  strength,  and  duetility  of  the  fibrils 
make  the  dominant  eontribution  to  the  eorresponding  fiber’s 
and,  in  turn,  yam’s  properties.  As  far  as  the  transverse  strength 
of  fibers  and  yams  are  eoneemed,  they  are  expeeted  to  be 
eontrolled  by  weak  interfibril  (and  in  the  ease  of  yams,  weak 
interfiber)  van  der  Waals/meehanieal  interloeking  forees  and 
to  be  less  dependent  on  the  fibril’s  transverse  meehanieal 
properties. 

The  results  obtained  in  Ref  10  ean  be  summarized  as 
follows:  (a)  fibril’s  stiffness  was  found  to  be  affeeted  by  both 
the  type  and  eoneentration  of  the  defeets.  An  example  of  the 
typieal  results  pertaining  to  the  effeet  of  defeet  type/eoneen- 
tration  on  the  Kevlar  fibril’s  stifftiess  is  depieted  in  Fig.  4(a) 
and  (b).  In  Fig.  4(a),  the  ease  of  -COOH  ehain  end  defeet  is 
eonsidered,  while,  in  Fig.  4(b),  the  ease  of  -SO3H  side  group 
defeet  is  analyzed;  (b)  As  far  as  the  fibril’s  axial  strength  is 
eoneemed,  it  was  found  to  be  affeeted  by  defeet  type  and  the 
size  and  eomposition  of  the  largest  defeet  eluster.  Examples  of 
typieal  results  obtained  are  shown  in  Tables  2  and  3.  The 
elusters  eomposed  of  defeets  of  the  same  type  are  depieted  in 
Table  2,  while  the  effeet  of  mixed-type  defeet  (ehain  end  and 
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defect 


Fig.  2  Examples  of  the  molecular-level  defects  in  Kevlar®  fibers:  (a)  a  COOH  chain-end  defect;  (b)  a  NH2  chain-end  defect;  (c)  a  COOH 
chain-end  defect  cluster;  (d)  a  SO3H  side-group  defect;  (e)  a  SO3H  side-group  defect  cluster;  (f)  a  microvoid  defect;  (g)  a  H2SO4  interstitial 
defect;  and  (h)  a  collection  of  chain  ends  forming  a  defect  band 


side  group)  clusters  is  depicted  in  Table  3.  An  examination  of 
the  results  reported  in  Tables  2  and  3  reveals  that  the  extent  of 
the  loss  of  fibril’s  axial  strength  increases  with  the  cluster  size/ 
strength;  (c)  The  aforementioned  finding  that  the  fibril’s  axial 
strength  is  a  function  of  the  largest  defect  cluster  and  its 
composition  indicated  that  fibril’s  strength  is  a  stochastic 
property.  Hence,  in  Ref  10,  results  like  those  displayed  in 
Tables  2  and  3  are  combined  with  the  prototypical  defect 
concentration  densities.  Table  1,  to  derive  the  associated 
Kevlar®  fibril-strength  distribution  function.  Fig.  5(a)  and  (b); 
(d)  In  accordance  with  Fig.  3(a),  Kevlar®  fibers  are  assumed  to 


be  transversely  isotropic  (with  the  fiber  axis  being  the  unique 
material  direction)  although  the  constituent  fibrils  possess 
orthotropic  symmetry.  Consequently,  (i)  fiber  axial  stiffness  and 
strength,  and  Gi  are  set  equal  to  their  fibril’s  counterparts; 
(ii)  the  fiber  transverse  normal  stiffness,  E22  =  ^33,  is  set  equal 
to  the  more  compliant  E33  fibril’s  modulus;  (iii)  all  the  three 
shear  moduli  G12,  G13,  and  G23  are  set  equal  to  the  same  value 
which  is  controlled  by  the  compliant  Kevlar®  intersheet  sliding 
resistance  (i.e.,  G12  of  the  fibrils);  and  (e)  To  assess  the 
yam-level  stiffness  and  strength  properties  from  the  known 
fiber  properties,  additional  effects  arising  from  the 
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defect 


Fig.  2  continued 

fiber-twist-induced  interfiber  friction  had  to  be  taken  into 
account.  This  procedure  mainly  affected  the  yam  axial  strength 
(or  more  precisely,  yam  axial  strength  distribution).  It  was 
generally  found  that  the  fibril-level  axial  strength  distribution 
function  displayed  in  Fig.  5(a)  and  (b)  was  modified  in  the  yam 
case  by  an  interfiber  friction  dependent  (positive)  strength  shift. 
For  example,  for  a  prototypical  extent  of  fiber  twist,  an  increase 
in  the  strength  of  ca.  0.2  GPa  was  found. 

4.  Enrichment  of  Yen  (Ref  8)  Composite 
Materiai  Modei 

In  this  section,  some  of  the  findings/results  reported  in 
section  3  are  used  to  enrich  the  continuum-level  material  model 
developed  by  Yen  (Ref  8).  Specifically,  Kevlar®-fiber’s  strength 
distribution  functions  like  the  one  displayed  in  Fig.  5(a)  and  (b) 
are  used  to  impart  a  stochastic  character  to  the  Kevlar®-fiber’s 
tensile  strength  and  transverse- shear  strength  parameters 
SyT,  and  Sy^s-  It  should  be  recalled  that  these  strength 
parameters  are  treated,  both  in  the  material  model  developed  by 
Yen  (Ref  8)  and  in  its  present  modification,  as  non-degradable 
material  properties.  To  demonstrate  the  effect  of  material  model 


Defect  band 


modifications  on  the  material  mechanical  response,  the  case  of 
uniaxial  tension  in  the  weft  direction  is  first  considered.  This  is 
followed  by  a  dynamic  computational  analysis  in  which  a  rigid 
projectile  of  a  right  circular  cylindrical  geometry  is  impacting, 
at  a  zero  obliquity  angle,  a  target  laminate  of  rectangular- 
parallelepiped  geometry.  In  the  latter  case,  it  was  critical  to 
demonstrate  how  much  the  proposed  material-model  modifica¬ 
tions  can  affect  the  composite  laminate  computed  ballistic  limit, 
as  quantified  by  the  V50  (the  projectile  incident  velocity  at 
which  the  probability  of  target  penetration  is  50%). 

4.1  Uniaxial  Tension  Test 

In  this  section,  a  quasi-static  uniaxial  tension  test  is  simulated 
using  the  conventional  finite  element  method.  Details  of  the 
finite  element  procedure  used  are  identical  to  the  one  reported  in 
our  prior  study  (Ref  11,  20-22).  In  short,  a  rectangular- 
parallelepiped-shaped  composite-material  test  sample  with  its 
longitudinal  axis  aligned  in  the  weft  direction  is  subjected  to  a 
progressively  larger  axial  strain  (while  no  lateral  constraints  are 
imposed,  to  attain  a  uniaxial  stress  state  condition).  To  identify 
differences  in  the  tensile-response  of  the  test  sample  brought 
about  by  the  proposed  modifications  in  the  material  model, 
various  global-response  and  spatial  distributions  of  material 
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Fig.  3  Topological  relationships  between  (a)  fibrils  and  fibers;  and 
(b)  fibers  and  yams  adopted  in  the  present  investigation 


state  variables  are  monitored/generated.  Details  of  the  results 
obtained  are  presented  below. 

A  eomparison  of  the  uniaxial  tension  (loading)  true  stress 
versus  uniaxial  tension  true  strain  eurves  for  a  prototypieal 
Kevlar®  plain  woven  eomposite  material  obtained  using  the 
original  Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in 
the  present  study  is  shown  in  Fig.  6.  Examination  of  this  figure 
reveals  that 

(a)  almost  up  to  the  onset  of  failure,  the  two  eurves  are 
effeetively  eoineident;  (b)  However,  the  eomposite  updated 
material  model  shows  that  the  onset  of  failure  oeeurs  at  a  3-4% 
lower  stress  level  and  a  slightly  smaller  strain;  and  (e)  The 
updated  material  model  prediets  a  higher  rate  of  material 
damage/degradation.  In  other  words,  the  material  displays 
almost  a  brittle-like  behavior. 

In  Fig.  7(a)-(d)  and  8(a)-(d),  a  eomparison  is  made  between 
temporal  evolutions  of  spatial  distribution  of  the  axial  tensile 
stress  in  the  weft  direetion  in  the  same  prototypieal  Kevlar® 
plain  woven  eomposite  material,  obtained  using  the  original 
Yen  (Ref  8)  material  model  and  the  updated  material  model 
proposed  in  the  present  study.  An  examination  of  these  results 
reveals  that,  at  eaeh  level  of  the  imposed  overall  axial  tensile 
strain,  the  spatial  distribution  of  the  axial  tensile  stress  is  fairly 
uniform.  This  finding  is  somewhat  infiueneed  by  the  ehoiee  of 
the  stress  seale  (100-600  MPa)  used  in  these  figures.  In  other 
words,  to  obtain  distinet  eolors  for  speeifie  stress  levels  over  the 
entire  loading  eyele,  a  eonstant/eommon  stress  range  was  used. 

In  Fig.  9(a)  and  (b),  a  eomparison  is  made  between  spatial 
distribution  of  the  axial  tensile  stress  in  the  weft  direetion  in  the 
same  prototypieal  Kevlar®  plain  woven  eomposite  material, 
obtained  using  the  original  Yen  (Ref  8)  material  model  and  the 
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Fig.  4  The  effect  of  the  concentration  of  the  (a)  -COOH  chain  end 
group;  and  (b)  -SO3H  side  group  on  the  orthotropic  elastic  stiffiiess 
constants  of  the  Kevlar®  fibrils 


updated  material  model  proposed  in  the  present  study  at  the 
instant  when  the  overall/average  axial  stress  experienees  a 
maximum.  An  examination  of  the  results  displayed  in  these 
figures  shows  that,  while  in  the  ease  of  the  original  Yen  (Ref  8) 
material  model.  Fig.  9(a),  spatial  distribution  of  the  axial  stress 
is  indeed  quite  uniform,  in  the  ease  of  the  upgraded  eomposite 
material  model  the  axial-stress  distribution  is  somewhat 
nonuniform.  To  explain  these  differenees,  a  eomparison  is 
made  in  Fig.  10(a)  and  (b)  between  the  spatial  distributions  of 
the  extent  of  axial  tension/trans verse  shear  damage  in  the  weft 
direetion  in  Kevlar®  plain  woven  eomposite  material  using  the 
original  Yen  (Ref  8)  material  model  and  the  updated  material 
model  proposed  in  the  present  study  at  the  instant  when  the 
overall/average  axial  stress  experienees  a  maximum.  It  is  seen 
that,  while  in  the  ease  of  original  Yen  (Ref  8)  material  model 
this  damage  variable  is  fairly  uniformly  distributed  (and  equal 
to  zero),  in  the  ease  of  the  updated  material  model,  one  of  the 
ten  eubie  finite  elements  eonstituting  the  tensile  speeimen  has 
aequired  a  non-zero  value  of  this  damage  variable.  This  finding 
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Table  2  A  summary  of  the  effect  of  the  most  common  defects  found  in  Kevlar®-based  materials  and  their  (same  defect- 
type)  clusters  on  the  PPTA  fibril’s  axial  tensile  strength  (in  GPa) 


Flaw  group 

Flaw 

Number  of  defects  in  the  largest  cluster 

1  2  3 

Isolated  ehain  ends  (point  defeet) 

-COOH 

6.0 

5.6 

5.0 

-COO“Na+ 

-NH2 

6.1 

5.8 

5.2 

-NH3^HS04“ 

Side  groups  (point  defeet) 

-SO3H 

5.8 

5.1 

4.1 

-S03“Na+ 

Voids  and  interstitials  (point  defeets) 

Mierovoids 

3.8 

N/A 

N/A 

Mobile  Trapped  H2SO4 

3.9 

N/A 

N/A 

Defeet  bands  (planar  defeets) 

NH3^HS04“  agglomerated  ehain  ends 

3.2 

N/A 

N/A 

Table  3  A  summary  of  the  effect  of  the  composition  and 
the  size  of  the  mixed  S03H/S03“Na^based  side  group  and 
COOH/COO“Na^  based  chain  end  defect  clusters  on  the 
PPTA  fibril’s  axial  tensile  strength  (in  GPa) 


Number  of  SO3H/SO3  Na^  defects  in  the 
cluster 


1 

2 

3 

Number  of  COOH/COO“Na^ 

defeets  in  the  eluster 

1  5.5 

4.8 

3.8 

2  5.3 

4.4 

3.4 

3  4.5 

3.8 

2.7 

explains  not  only  the  axial-stress  nonuniformities  observed  in 
Fig.  9(b)  but  also  the  more  brittle  nature  of  the  eomposite 
material.  In  other  words,  onee  one  of  the  elements  aequires  a 
non-zero  value  of  the  uniaxial  tension/transverse  shear  damage 
variable,  further  deformation  beeomes  loealized  within  this 
element.  This  eauses  both  the  struetural  weakening  of  the 
material  in  this  element  (due  to  loss  of  the  assoeiated  Young’s 
modulus)  and  the  geometrieal  weakening  (due  to  a  reduetion  in 
the  eross-seetional  area  supporting  the  applied  axial  load).  It 
should  be  noted  that  the  observed  differenees  between  the 
results  presented  in  Fig.  9(a)  and  (b),  as  well  as  the  differenees 
between  the  results  presented  in  Fig.  10(a)  and  (b),  are  due  to 
the  stoehastie  nature  of  material  strength  parameters  imple¬ 
mented  into  the  updated  material  model.  In  other  words,  the 
failed  element  identified  in  Fig.  9(b)  and  10(b)  is  assoeiated 
with  the  lowest  level  of  tensile  strength. 

The  results  displayed  in  Fig.  6-10  elearly  reveal  that 
ineorporation  of  the  realistie  miero-meehanies  effeets  sueh  as 
the  operation  of  various  fiber/matrix-damage  meehanisms  and 
the  stoehastie  nature  of  the  material  strength  ean  signifieantly 
alter  the  global  as  well  as  the  loeal  response  of  the  material  to 
the  imposed  loading.  While  this  finding  ean  be  eonsidered  as 
expeeted,  its  direet  quantitative  experimental  validation  is  quite 
ehallenging  and  is  beyond  the  seope  of  the  present  study. 

4.2  Laminate  Transverse  Impact  by  a  Rigid  Cylindrical 
Projectile 

In  this  seetion,  the  problem  of  eomposite  laminate  transverse 
impaet  by  a  right  eireular  eylindrieal  rigid  projeetile  is  analyzed 


Fig.  5  (a)  Kevlar®-fibril’s  axial  strength  probability  density  fune- 

tion;  and  (b)  the  eorresponding  eumulative  distribution  funetion 
determined  in  the  present  study  using  moleeular-type  eomputational 
analyses 
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Fig.  6  A  comparison  of  the  uniaxial  tensile  true  stress  vs.  uniaxial 
tensile  true  strain  behavior  (in  the  weft  direction)  in  a  prototypical 
Kevlar®  plain  woven  composite  material  obtained  using  the  original 
Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in  the  present 
study 
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Fig.  8  Spatial  distribution  of  the  axial  tensile  stress  in  the  weft 
direction  in  a  prototypical  Kevlar®  plain  woven  composite  material 
obtained  using  the  upgraded  model  proposed  in  the  present  study  at 
the  average  axial  strain  levels  of  (a)  0.004;  (b)  0.014;  (c)  0.024;  and 
(d)  0.032 
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Fig.  7  Spatial  distribution  of  the  axial  tensile  stress  in  the  weft 
direction  in  a  prototypical  Kevlar®  plain  woven  composite  material 
obtained  using  the  original  Yen  (Ref  8)  material  model  at  the  aver¬ 
age  axial  strain  levels  of  (a)  0.004;  (b)  0.014;  (c)  0.024;  and  (d) 
0.033 

using  a  dynamic  explicit  finite  element  eomputational  analysis. 
Sinee  details  of  the  finite  element  proeedure  used  is  identieal  to 
the  one  reported  in  our  prior  study  (Ref  11),  they  will  not  be 
repeated  here. 

An  example  of  the  prototypieal  results  obtained  in  the 
present  study  is  given  in  Fig.  1  l(a)-(d).  The  results  displayed  in 
Fig.  ll(a)-(d)  pertain  to  the  spatial  distribution  of  the  eompos- 
ite-laminate  material  during  impaet.  It  should  be  noted  that  to 
enable  monitoring  of  the  penetration  proeess,  the  target  plate  is 
made  transparent.  The  results  displayed  in  Fig.  ll(a)-(d)  elearly 
reveal  the  development  of  the  baek  faee  bulge.  In  general,  it 
was  found  that  details  regarding  the  spatial  distribution  and  the 


589MPa 


(a) 


Fig.  9  A  comparison  between  the  spatial  distributions  of  the  axial 
tensile  stress  in  the  weft  direetion  in  Kevlar®  plain  woven  eomposite 
material  at  the  instant  when  the  overall/average  axial  stress  experi- 
enees  a  maximum,  obtained  using  (a)  the  original  Yen  (Ref  8)  mate¬ 
rial  model;  and  (b)  the  updated  material  model  proposed  in  the 
present  study 
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Fig.  10  A  eomparison  between  the  spatial  distributions  of  the  ex¬ 
tent  of  axial  tension/transverse  shear  damage  in  the  weft  direetion  in 
Kevlar®  plain  woven  eomposite  material  at  the  instant  when  the 
overall/average  axial  stress  experienees  a  maximum,  obtained  using 
(a)  the  original  Yen  (Ref  8)  material  model;  and  (b)  the  updated 
material  model  proposed  in  the  present  study 

temporal  evolution  of  the  target  material  during  impaet  are 
affeeted  by  the  modifieations  made  to  the  original  Yen  (Ref  8) 
eomposite  material  model.  Details  regarding  the  effeet  of  these 
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Fig.  11  Typical  results  pertaining  to  the  temporal  evolution  of  the 
spatial  distribution  of  the  composite-laminate  material  during  impact 
by  a  right-circular  cylindrical  rigid  projectile 
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Fig.  12  A  typical  projectile  residual  velocity  vs.  projectile  incident 
velocity  plot  showing  differences  in  F50  as  predicted  by  the  original 
Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in  the  present 
study 


material  model  modifieations  on  the  temporal  evolution  and 
spatial  distribution  of  the  target-material  and  its  various 
meehanieal  and  mierostruetural  fields  will  be  presented  in  our 
future  eommunieation.  In  the  remainder  of  this  seetion,  results 
pertaining  to  the  effeet  of  material  model  modifieations  on  V50 
are  presented  and  diseussed. 


Typieal  results  pertaining  to  the  effeet  of  projeetile  initial 
veloeity  on  the  projeetile  residual  veloeity  in  the  ease  of  the 
original  Yen  (Ref  8)  eomposite  material  model  and  its  present 
modifieation  are  displayed  in  Fig.  12.  It  is  seen  that  V50 
(defined  here  as  the  projeetile  initial  veloeity  at  whieh  the 
projeetile  residual  veloeity  is  zero)  predieted  by  the  analysis  in 
whieh  the  modified  eomposite  material  model  is  used  is  ea.  5% 
lower  than  its  eounterpart  assoeiated  with  the  use  of  the  original 
Yen  (Ref  8)  eomposite  material  model.  This  effeet  is  deemed  as 
being  signifieant  and  justifies  the  present  efforts  aimed  at 
ineorporating  fine-seale  mierostruetural  effeets  into  the  eontin- 
uum-level  material  model.  However,  no  relevant  experimental 
data  are  presently  available  to  be  able  to  judge  the  potential 
improvements  in  the  experiment/ealeulation  agreement  brought 
about  by  the  use  of  the  modified  eomposite  material  model. 
This  defieieney  in  the  relevant  experimental  data  will  be 
addressed  in  our  future  study. 


5.  Summary  and  Conclusions 

Based  on  the  results  obtained  in  the  present  study,  the 
following  summary  remarks  and  main  eonelusions  ean  be 
drawn: 

1 .  A  eontinuum-level  rate-dependent  material  model  for 
plain  woven  Kevlar®-fiber-reinforeed  polymer  matrix 
eomposite  has  been  upgraded  by  ineorporating  the  effeets 
of  fibril/fiber/yam  level  mierostrueture  phenomena  and 
proeesses. 

2.  Speeifie  mierostruetural  effeets  ineorporated  into  the  eon¬ 
tinuum-level  material  model  aeeount  for  the  role  of  vari¬ 
ous  material  synthesis/fiber  proeessing-indueed  defeets  on 
the  material  strength,  or  more  speeifieally  on  the  material 
strength  distribution  fimetion. 

3.  For  a  prototypieal  plain  woven  Kevlar®-fiber-reinforeed 
polymer  matrix  eomposite,  the  results  obtained  using  a 
series  of  finite  element  analyses  show  that  the  predieted 
response  of  the  material  both  under  quasi-statie  and 
dynamie  ballistie-impaet  eonditions  is  signiheantly  affeeted 
as  a  result  of  the  ineorporation  of  mierostruetural  effeets. 
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